Kohn-Luttinger instability of the t-t' Hubbard model in two dimensions: 

variational approach 
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An efTective Hamiltonian for the Kohn-Luttinger superconductor is constructed and solved in the 
BCS approximation. The method is appHed to the t-t' Hubbard model in two dimensions with the 
following results: (i) The superconducting phase diagram at half filling is shown to provide a weak- 
coupling analog of the recently proposed spin liquid state in the J1-J2 Heisenberg model, (ii) In the 
parameter region relevant for the cuprates we have found a nontrivial energy dependence of the gap 
function in the dominant d-wave pairing sector. The hot spot effect in the angular dependence of 
the superconducting gap is shown to be quite weak. 



I. INTRODUCTION 

In a seminal paper, Kohn and Luttinger0 (KL) pre- 
sented a perturbative argument showing that the generic 
ground state of a Fermi liquid at weak coupling is super- 
conducting, even if the bare electron-electron interaction 
is not attractive in the Cooper channel in any angular 
momentum sector. Recently, this result has been red- 
erived utilizins^ various versions of the renormalization 
group method.D'DlJ 

However, no attempt has been made so far to con- 
struct a variational wavefunction for a KL superconduc- 
tor, which could be (at least in principle) tested in a 
direct numerical simulation. The purpose of this pa- 
per is to construct such a wavefunction for the repulsive 
Hubbard model in the limit of weak coupling U ^ W, 
where U is the local repulsive interaction and W is the 
electron bandwidth. Our strategy follows quite closely 
the canonical transformation approach-|to the electron- 
phonon problem introduced by FrohlichB and sketched in 
the Appendix. In some sense it is also similar to the cus- 
tomary analysis of the Hubbard model close to half filling 
in the limit of strong coupling U ^ W. Namely, in 
latter case a canonical transformation can be founc 
which transforms the purely repulsive Hubbard model to 
the so-called t-J model which explicitly contains an at- 
tractive interaction (of the order J ^ W'^/U) favoring the 
formation of singlet superconductivity. Unfortunately, 
the t-J model contains a constraint on the local number 
of electrons and as such is very difficult to study. Here, 
on the other hand, we construct a canonical transforma- 
tion which eliminates the scattering of quasiparticles to 
first order in U/W and generates a weak-coupling model 
with an attractive interaction of the order U'^/W. 

The outline of the paper is as follows. The construction 
of the effective model is described in Section II and our 
method is applied to the two-dimensional t-t' Hubbard 
model in Section III. 



described by the Hamiltonian H = Hq -\- Hi-\- H2, where 
Hq = ^ Ek'T'k.o- is the kinetic energy operator and the 
interaction term has been split into two parts: 
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where N = N-^ -\- and S are the total electron num- 
ber and the total spin, respectively. Terms which are 
not extensive in the thermodynamic limit L — > 00 have 
been neglected. H2 is the generic interaction term which 
scatters electrons from ki t to ka t and from k2 i to 
k4 J,. The prime on the summation means that terms 
with ki = ks (forward scattering channel), ki = k4 (ex- 
change channel), and ki + k2 = (Cooper channel) are 
excluded from H2 and are singled out into Hi. The op- 
erator Hq -I- Hi can be thought of as a reduced BCS 
Hamiltonian of a Landau Fermi liquid, while H2 con- 
tains the scattering processes leading to a finite lifetime 
of the Landau quasiparticles. 

Now we look for a canonical transformation from the 
bare electrons to the dressed ones, H — e^^He~'^^ , such 
that the scattering of the quasiparticles vanishes to first 
order in U . This happens if H2 + i[S, Hq] = in which 
case we can write H = Hq -\- Hi -\- i[S, Hi] -\- i[S, i?2]/2, 
where terms of order 0{U^) have been neglected. Making 
use of the identity 
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where the indices P, Q, i?, S include also the spin label, 
one verifies readily that 
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II. EFFECTIVE HAMILTONIAN 

We consider the Hubbard model on a square lattice 
with periodic boundary conditions and L ~ I x I sites. 



is the Hermitian generator we are looking for. Note that 
the forward, exchange, and Cooper channel processes do 
not contribute to S. In order to calculate the commuta- 
tors [S, Hi] and [S, H2] we first note that the operators 
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A = C3|Ci| and B = cJ^|Cq| commute with the operators Combining Eqs. (^|j^) we thus find the identity 
X = c^4^^C2i and Y ~ Cg^cp^ and therefore 

[AX, BY] = [A, B]XY + BA[X, Y]. (5) 



making use of which one can straightforwardly calculate the effective Hamiltonian H . The result is fairly involved, 
including also three-body forces among the electrons. Nevertheless, the expectation value E = {tp\H\tp) in a BCS 
state ji/j) — nk(uk + ''^kC|{(.|C^i(.|)|0) takes a simple and physically transparent form, 

E = 2^ Ek/k. + J^ + j2^ Vkpb^bp + j22^ _ 4i+k.,k3+k4, (6) 

ka ^ * ^ k,p ^ {k} + ^''^ '^'^^ 



where we have introduced a momentum-resolved BCS 
order parameter bp — {^^:\c-plCp^\^p) and we assumed 
{tfj\S'^\tp) = 0. In the Cooper channel the effective in- 
teraction reads Vkp = U + C/^Xi(k -f p, £p — £k), where 
Xi(q, uj) is the real part of the particle- hole susceptibility 
X(q,t^) = L-^ Ek(/k - /K+q)/(eK+q - £k - - iO). 
The physical content of the KL argument is that the ef- 
fective interaction between the dressed electrons can be 
attractive in the Cooper channel in some angular mo- 
mentum sector. If this is the case, then the energy 
E — {iple^^' He~''^ \ip) is minimized by nonvanishing bp 
and the wavefunction [ip) — e~''^|'(/') provides a varia- 
tional ansatz for the KL superconductor. 

Before proceeding we would like to point out that our 
variational approach bears some similarity to. a very re- 
cent formulation of the flow equation method.B The com- 
parison of both methods is left for future. 



III. BCS VARIATIONAL SOLUTION OF THE 
EFFECTIVE HAMILTONIAN 

We wish to apply the above formalism to the study 
of the t-t' Hubbard model in two dimensions whose sin- 
gle particle dispersion is Ek — ~2t{coskx + cos ky) + 
At' cos cos ky , where we have set the lattice constant 
a — 1. The superconducting phase diagram of the weak- 
coupling t-t' Hubbard model has been studied previ- 
ously by considering the effective KL interaction — 
U + U^^i (k -I- p, 0) with both k and p lying at the Fermi 
surface.El Roughly speaking, such an approach provides 
us with the coupling constant g, but not with the prefac- 
tor in the BCS formula for the transition temperature, 
T, = Oexp(-l/g). 

In what follows we perform the standard BCS mini- 
mization of an effective free energy whose T = limit re- 
duces to the variational energy Eq. M). We assume that 



the interaction U is sufficiently weak so that no particle- 
hole instabilities can develop. (In particular, this means 
that we can't study systems close to the van Hove fill- 
ing.) Concentrating on lowest order effects in U/W, we 
neglect the last term in Eq. (^) which represents the stan- 
dard second-order perturbation theory effects, and treat 
it as an additive constant. A more detailed study of 
Eq. taking into account also the last term is post- 
poned to future publications. Moreover, since we expect 
that Tc <SC t, we make use of the T = Fermi functions in 
the calculation of the particle- hole susceptibility x(q, t^)- 
A standard BCS calculation leads to the self-consistent 
equation for the gap function Ak = J^p ^kp^p, 

P P 

where Ep = {^^ + Ap)^/^, ep — ^, and fj. is the 

chemical potential. 

For temperatures infinitesimally below the mean field 
transition temperature, the gap equation can be lin- 
earized with respect to Ap and can be written as L'k = 
Ep-fkp(T)£)p, where Dy, = Ak^kT-^/^ ^ dimen- 

sionless gap function, 0k = [tanh(^k/2T)/2^k]^^^, and 
/kp(T) = — i^^Vkp'/'k^p is a real symmetric matrix cor- 
responding to a dimensionless pair scattering function. 
At high temperatures /kp(T) oc and the gap equa- 
tion does not have solutions. With decreasing T the 
maximal eigenvalue A of /kp (T) grows and the mean field 
transition temperature can be calculated from the con- 
dition A(Tc) 1. 

We have studied the eigenvalues of Ikp{T) numeri- 
cally. In order to minimize the cost of the calculation, 
the matrix /kp was calculated on a special lattice, see 
Fig. |l]. The lattice consists of a sequence i — 1,2, ... ,n 
of n regular li x li grids, which become progressively 
more dense as the distance to the Fermi level dimin- 
ishes, Zi+i = 2li. The borders between the subsequent 
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FIG. 1: An example of the special set of k points in the first 
Brillouin zone for t' /t = 0.3 and p — 0.8. Under approaching 
the Fermi energy, the grid changes in each level (there are 4 
levels in total) from 32 x 32 to 256 x 256. 



levels are constructed as follows: Let the maximum 
and minimum band energies be ^max and ^min and let 
^0 = niax(^max, I'fminD- Wc construct a sequence of en- 
ergies = g'^o with a quotient q < I. For i < n, the 
i-th level grid li x 1^ is realized for all k points satisfy- 
ing < l^kl < Ci-ii whereas the n-th grid applies for 
ICkl < £,n~i- The lattice is characterized by three param- 
eters: q, h, and the number of levels n. In what follows 
we always take q = and li = 32, and instead of the 
number of levels we describe the lattice by specifying the 
finest k-point grid Z„ x 

The susceptibility x(q, was calculated on usual lat- 
tices X /„ by a straightforward modification of the 
method described in Ref. |l^ which makes use of the 
fast Fourier transform algorithm. We have used 4Z„ time 
points and the Nyquist frequency was chosen to be 814^. 

The matri*,Jkp was diagonalized by the modified Lanc- 
zos method.tll The initial vector was generated randomly 
in one octant (quadrant) of the Brillouin zone and sub- 
sequently continued to the whole zone so as to transform 
according to the one-dimensional s, d, d^y and g (two- 
dimensional p) irrediLcible representations of the point 
group of the square. E£l 

As an illustrative example, the temperature depen- 
dence of the maximal eigenvalue in the dxy sector for 
t' /t = 0.78 and electron density p = 1 is shown in Fig. |^ 
for various In x Z„. Several points are worth mention- 
ing. For a finite lattice size In, the divergence of the 
particle-particle susceptibility is cut off at temperatures 
T* ~ vp/ln, and therefore A(T) saturates at low temper- 
atures T <C T*. Thus, even making use of a moderate 
interaction strengthllj U — W/2 (this value of U is used 
in most numerical examples) the largest eigenvalue A of 
Jkp is typically less than unity down to the lowest tem- 
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FIG. 2: The temperature dependence of the largest eigenvalue 
in the d^jz-wave sector for U = W/2, t' /t = 0.78, and p = 1. 
The line is a fit to A(r) = glog{Q/T). The fitting parameters 
g = 0.04 and n/t = 480 imply Tc/t = 6.05 x 10"^ 



peratures. Therefore we can't determine Tc directly from 
A(Tc) — 1. Instead, the T-dependence of the eigenvalue is 
fitted in the low temperature limit by a BCS-like expres- 
sion A(T) ~ g\og{n/T) and afterwards Tc is calculated 
from Tc — Qexp(—l/g). This procedure is shown explic- 
itly in Fig. ||. It should be pointed out that usually we 
can determine the transition temperature with a reason- 
able degree of confidence only if Tc/t > 10^^. In this 
sense the parameter values used in Fig. are somewhat 
special. 

Let us also stress that the large value of estimated 
from Fig. || should not be interpreted as a large pair- 
ing energy scale. This is evident from an equivalent but 
physically more meaningful expression A(T) = A(To) + 
glog{To/T) which makes explicit use of the pairing en- 
ergy To (its value is arbitrary, but can be estimated from 
the maximal temperature where the logarithmic scaling 
of A applies, Tq ^ 0.005t in our example). The large 
value of = Tq cxp[X{To) / g] is a consequence of the 
large A(ro)j which in turn measures the contribution of 
the high energy processes to A. 

A. Half-filled t-t' Hubbard model 

Now let us apply the abovementioned method to the 
study of the superconducting phase diagram of the t-t' 
Hubbard model at half filling {p — 1). In the strong- 
coupling limit U ^ W, this would correspond to the 
Hciscnberg model H ^ Ji Y,{i.j) • Sj -t- J2 ^i-Sj, 
where (i, j) and ((i, j)) are pairs of nearest neighbor and 
next-nearest neighbor sites, respectively, Ji = At'^/U, 
and J2 = 4(t')^/C/. There is general consensus that, 
due to frustration, for 0.4 < J2/J1 < 0.6 the Heisenberg 
model does not exhibit the Neel order (see, e.g. Ref. [l4| 
and references therein). The nature of the emergent state 
has been discussed intensively, but very recently it has 
been shownllil that a projected BCS wavefunction with a 
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FIG. 3: Superconducting phase diagram for the half filled 
Hubbard model at [/ = W 12. The inset shows the coupling 
constant g in the d and s sectors, calculated according to 
Ref. 1^ in the region where Tc is not directly accessible. 



Then, if we denote that value of g2 for which the cross- 
ing happens as 17201 it can be shown readily that a mixed 
d + idxy pairing state is in fact stabilized at zero tem- 
perature in a finite window of the control parameter for 
which 5^0^ - 1/2 < g^^ < g^Q + 1/2. Note that in the 
weak-coupling limit this window is quite narrow. 

Before proceeding we should mention that the con- 
cept of mixed-symmetry superconductivity has been in- 
troduced previously in the context p£ strong coupling 
models, for the t-J model by KotliaiO and £of the t-t'- 
J1-J2 model in a recent paper by Sachdev.Ej Our dis- 
cussion shows that (as pointed out already in RefJ|), 
mixed-symmetry superconductivity is in fact a generic 
consequence of phase diagrams like the one shown in 
Fig. y, where superconducting states of different sym- 
metry cross. 



B. d-wave region of the t-t' Hubbard model 



mixed d and d^y symmetry provides a very good approx- 
imation for the J1-J2 model with J2/J1 close to 1/2. 

In this paper we look for a weak-coupiing analogue of 
the phase proposed by Capriotti et alS-i To this end, in 
Fig. y we plot Tc in the d-wave and d^y-wave sectors as a 
function of a = 2t'/{t + 2t') for U = W/2 and p = 1. The 
dominant pairing instability is seen to change from the d- 
wave at small t' /t to the dxy symmetry at large t' /t, since 
the susceptibility x(q, 0) peaks at (tt, tt) and (7r,0) in 
these two limits. Surprisingly, in the region correspond- 
ing to J2/ Ji « 1/2, namely t' /t « 2'^^'^ or a « 0.59, we 
find that neither d nor d^^y-wave pairing instabilities are 
dominant, and a small island of s-wave pairing is realized. 
For a E [0.45, 0.57] the transition temperatures are too 
low to be accessible by the present method. Therefore 
we show in Fig. ^ the coupling constant g in the (domi- 
nant) d and s-wave symmetry sectors calculated accord- 
ing to Ref. ^. From the qualitative agreement of these 
two different methods we conclude that in the vicinity of 
a ^ 0.55 the pairing symmetry changes from d to s-wave. 

We expect that for a ~ 0.59, upon increase of the inter- 
action strength there is presumably a true Mott-Hubbard 
superconductor-insulator transition at some critical value 
of U. We hypothesize that before this happens, with in- 
creasing U the d and d^y phases grow at the expense of 
the s phase, until at some stage the s phase disappears 
from the phase diagram. On the other hand, underneath 
the crossing point of the d-wave and d^y-wave Tc's, we 
expect a region with mixed d + id^y pairing. Such mixed- 
symmetry pairing can be justified by a simple model cal- 
culation for a two dimensional isotropic system with a 
separable interaction V{(j),(j>') — Vi{(f>,(j)') + V2((/), </>') in 
the Cooper channel. Here Vi((/), 0') oc cos 20 cos 20' with 
coupling constant gi and cutoff wi leads to d-wave pair- 
ing and 1/2(0,0') oc sin 20 sin 20' with coupling constant 
52 and cutoff 0^2 yields d^y-yfa^ve pairing. For simplicity 
we further assume that by changing a control parameter 
{f /t in our case), only the coupling constant g2 changes. 



As another application of the present formalism, in the 
rest of this paper we concentrate on the d-wave region 
of the supexconducting phase diagram of the t-t' Hub- 
bard modelEI which is relevant to the cuprates. Close to 
the van Hove deasity, superconductivity competes with 
the SDW statectBEa and the question whether one of 
these two states or even more exotic statesEII wins is 
still open. Here we focus on a different question, namely 
that of the full k-dependence of the superconducting gap 
Ak away from the van Hove density. Fig. || shows the 
self-consistent solution Ak of Eq. (|^) for t' /t — 0.3, 
p — 0.8, and T = 0. The maximal gap on the Fermi 
surface A^ax = 0.0023t and the critical temperature 
T, = O.OOlSi imply A^^^JT^ « 1.8. 

The overall shape of Ak is well described by the sim- 
plest formula consistent with d-wave symmetry, Ak oc 
coskx — cosfcy. Nevertheless, fine structure on top of 
this overall shape is clearly visible. This is shown more 
quantitatively in Fig. ||, where Ak is plotted along two 
ky =const cuts of the Brillouin zone. A nontrivial struc- 
ture is seen to develop in the direction perpendicular 
to the Fermi surface which is analogous to the strong- 
coupling effects in conventional superconductors. 

It is well known that in the case of phonon mediated 
superconductivity the energy dependence of Ak leads to 
a structure in the density of states. In the Appendix we 
have shown that our method is capable to describe such 
effects. Therefore we have looked for features in the den- 
sity of states caused by the structure in Ak in our purely 
electronic model as well. The electron spectral function 
A(k, oj) and the density of states N{uj) have been calcu- 
lated from 

A{k,u;) = ul5{Lu~E^) + vlS{io + E^), (8) 
Niiv) = L-i^A(k,c.), (9) 

k 

where u^., v'^ = 2~^(l±^k/£'k)- Fig. ^ shows that at weak 
coupling (e.g. for U = W/2 when A,Tc <C t) no struc- 
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FIG. 4: The self-consistent solution Ak in the d-wave sector 
(in units of t) as a function of momentum in the first Brillouin 
zone for r = 0, (7 = W/2, t' /t = 0.3, and p = 0.8. Calculated 
on a full lattice L — 256 x 256. The momenta kx and ky are 
measured in units of tt. 




FIG. 5: Cuts of the data in Fig. ^ along the lines ky — tv (top) 
and ky = 5tt/8 (bottom). The crossing points of the ky — 
const lines with the Fermi surface are indicated by arrows. 
The inset shows the Fermi surface for t' /t — 0.3 and p = 0.8, 
and the location of the cuts. 



tures associated with the energy dependence of can be 
observed. This is because features in the energy depen- 
dence of Ak can be visible also in the density of states 
only if they appear in a k-point for which the single- 
particle energy is not much larger than Ak. In other 
words, in order to be observable in iV(a;), the structure 
has to appear at a distance of at most Sk ~ A/vp from 
the Fermi line. At weak coupling this criterion is not 
satisfied and therefore Fig. ^ does not exhibit any devia- 
tions from the textbook form of the density of states for 
a d-wave superconductor. 

We have solved the self-consistent equation Eq. (0) also 
for a larger interaction strength U = 0.75Vl^, when the 
gap increases by more than an order of magnitude with 
respect to U = W/2. In this case we did find a nontrivial 
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FIG. 6: Tunneling density of states for U = W/2, t' /t = 0.3, 
and p = 0.8. The superconducting gap obtained on a full 
lattice L = 256 x 256 was linearly interpolated to a lattice 
L = 8192 X 8192, over which the k-summation in Eq. ^ is 
taken. The delta functions appearing in the electron spectral 
function Eq. ^ were given a finite width 7/t = 1 x lO"*. 
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FIG. 7: Tunneling density of states for U = 0.75W, t' /t = 0.3, 
and p = 0.8. Calculated by the same method as in Fig. ^ with 
7/t = 5 X 10-''. 



structure in the density of states (see Fig. [t]), which could 
however be attributed to the van Hove singularity in the 
density of states of the noninteracting spectrum rather 
than to a many body effect. 

Now let us turn to the angular dependence of Ak- In 
Fig. ^ we plot A(i^)/A(0) along the Fermi line for T — 
and T/t = 1. The variable (p is the angle between the ra- 
dius vector of the Fermi surface point under study [mea- 
sured with respect to M = (tt, tt)] and the x axis. The 
overall shape of the high temperature data is close to the 
simple cos kx — cos ky form of the gap. With decreas- 
ing temperature Ak is seen to be locally more and more 
enhanced close to the crossing point of the Fermi line 
with the magnetic zone (hot spot). Nevertheless, this 
so-called hot spot effect is quite weak down to T — 0. 
Fig. ^ shows explicitly that the hot spot effect-is in qual- 
itative agreement with an earlier formulationLi in which 
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FIG. 8; Angular variation of a normalized superconducting 
gap along the Fermi line. Calculated on a full lattice L = 
256 X 256 for U = W/2, t' /t = 0.3, and p = 0.8. Top to 
bottom lines at if — 30°: T/t = 1, the simplest d-wave model 
function cos fc^; — cos ky, and T = 0; the lowest curve is a result 
obtained within the approach of Ref. H. 



ing state at t' /t ~ 0.7 is a weak-coupling analog of the 
wavefwuiction proposed for the J1-J2 model by Capriotti 
et alM 

As another applination, in the d-wave region of the 
t-t' Hubbard modeltl we have found a nontrivial en- 
ergy dependence of the gap function reminiscent of 
strong-coupling effects in conventional superconductors. 
There is a hot spot effect in the angular dependence of 
Ak, which is quite consistent with the estimates from 
calculationda in which the gap lives only on the Fermi 
line. 
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APPENDIX 



the energy dependence of the Cooper channel interaction 
is neglected and Ak lives cpiy on the Fermi line. On the 
other hand, a recent papeillH reports (within the formal- 
ism of Ref. P) a much stronger hot spot effect than found 
here. 



IV. CONCLUSIONS 

In conclusion, we have constructed a canonical trans- 
formation of the Hubbard model which eliminates quasi- 
particle scattering processes to first order in the inter- 
action U. The resulting effective Hamiltonian contains 
terms of order U'^/W where W is the bandwidth, which 
are attractive in the Cooper channel. This allows us 
to construct a variational wavefunction for the Kohn- 
Luttinger superconductor and to show explicitly that 
even a purely repulsive model gains energy by developing 
superconducting correlations. 

As an application of the method, the effective Hamil- 
tonian has been solved in the BCS approximation and 
the superconducting phase diagram of the t-t' Hubbard 
model at half filling has been found. The superconduct- 



As a test of the applicability of the variational method 
to an effective electronic Hamiltonian, we apply it to the 
case of phonon mediated superconductivity. We consider 
the simplest model of electrons with a local coupling to 
a phonon mode. The Hamiltonian of the system can be 
written H = Hq + H' , where 

Ho = X!^kCk,£TCk,fT + X!^q"q"q' 

k,(T q^O 

^' = ^ E 4+q..Ck,.k + aO. (A.l) 

V-^k,<T,q#0 

Following Frohlichji we seek a canonical transformation 
H — e'-^ He~^^ which eliminates electron-phonon cou- 
pling to first order in D. This is accomplished by 

/ 4+q,^Ck,^aq ^ 4+q,^Ck,^aLq \ 
^kT^^O ek+q - £k + t^q y ' 

(A.2) 

since, as one verifies easily, [iS,HQ\ = —H'. The effective 
Hamiltonian can therefore be written as H = Hq -\- H' 
where, to second order in D, 



H'^-[^S,H'] = 



k,k^.wto(-'^+<>--'^)^— ^ """" 



E E E^'^k-l-q.o-Ck-^q'.o- cj^+q+q'^ffCk.o-) 



q^^Oq'^O k,cr 



Ck',(T'Ck_f,q o-Ck^CT 

aq(aq' +aLqO ^ aLq/(aq' +QLqO 

Ek+q — £k — Wq £k-|-q — £k — Wq 
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The expectation value of H' in the product of a BCS state with a phonon state diagonahzing Hq, jV'), simphfies to 



k,k' k,k ,o- q7=0 




FIG. 9: Cut of the superconducting gap for phonon mediated 
superconductivity along the line ky — tt. Calculated on a 
full lattice L = 256 x 256 for t' /t = 0.45, and p = 1.0. The 
phonon parameters are uJo/t ~ 0.2 and T/t = 0.02 and the 
electron-phonon coupling D/t = 0.51. The dashed vertical 
lines correspond to quasiparticle energies = and ±ljo- 
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FIG. 10: Density of states in the superconducting state N{uj) 
normalized by its value in the normal state Nn{tjj). Calcu- 
lated from the data in Fig. ^ making use of the same method 
as in Fig. pi The dashed line is the BCS approximation 
N{Lo)/Nn{LJ) = Lj/Vcj'^ - A2 with A = 0.059t. 



where n{uj) is the Bose distribution function. The first 
term describes the effective phonon mediated electron- 
electron interaction in the Cooper channel with Vk^k' = 
2£'^a;k-k'/[(^k — ^k')^ ~ '^k-k']' ^he second term is the 
standard second-order perturbation theory correction to 
the ground state energy of the electron system, and the 
last term with Sajq = ~2D^xi(c1j "^q) is the phonon en- 
ergy renormalization. 

In order to facilitate a numerical solution of the gap 
equation Eq. (0) for phonon mediated pairing, we have 
made use of a Cooper channel interaction 



Ek - Ek' + o-(^^k-k' + «r) ' 



regularized by a finite phonon lifetime F. The BCS gap 
equation has been solved numerically on the same t-t' 
square lattice as for the Hubbard model. The supercon- 
ducting gap obtained numerically for an Einstein mode 
with frequency ujo/t — 0.2 and T/t = 0.02 for electronic 
parameters t' /t = 0.45 and p = 1.0 and electron- phonon 
coupling D/t = 0.51 is shown in Fig. ^. Note that the 
gap changes sign in the vicinity of the phonon energy 
luq. This nontrivial energy dependence of Ak translates 
into a feature of the tunneling density of states in the 
same energy region, see Fig. |lO|. It should be pointed 
out that although both of the above results axa in qual- 
itative agreement with the Eliashberg theoryjHj there is 
one important difference of the present approach with re- 
spect to Ref. |o[ Namely, within the Eliashberg theory 
the features of A(cl') and N{uj) occur at energy u!o + A{Q) 
and not loq as in our case. This difference is presumably 
due to the fact that within the canonical transformation 
approach, it is the bare elect ron e nergies which enter the 
energy denominators in Eq. ( [A.2| ) and not the renormal- 
ized energies as in the Eliashberg theory. However, in a 
truly weak coupling theory this difference is marginal. 
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